Pair correlations in sandpile model: a check of logarithmic conformal field theory 
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We compute the correlations of two height variables in the two-dimensional Abelian sandpile 
model. We extend the known result for two minimal heights to the case when one of the heights 
is bigger than one. We find that the most dominant correlation log r/r* exactly fits the prediction 
obtained within the logarithmic conformal approach. 
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Conformal field theory has proved to be extraordinar- 
ily powerful in the description of universality classes of 
equilibrium critical models in two dimensions [1]. Crit- 
ical exponents, correlation functions, finite-size scaling, 
perturbations and boundary conditions, among others, 
have all been studied within the conformal approach, and 
thoroughly (and successfully) compared with numerical 
data. 

More recently, increased interest has been turned to- 
ward logarithmic conformal theories, as a larger class of 
conformal theories, interesting in its own right, but also 
as a description of certain non-equilibrium lattice mod- 
els. In particular, dense polymers [2, 3], sandpile models 
[4, 5] and percolation [3, 6] are lattice realizations of log- 
arithmic conformal theories. An infinite series of such 
lattice models have been defined in [3]. 

The logarithmic theories are however much lesser un- 
derstood than the more usual, non- logarithmic ones. 
This is due to their higher level of complexity, which 
somehow refiect the complexity of the associated lattice 
models. Indeed the models mentioned above all have in- 
trinsic non-local features. In this respect, it may appear 
to be strange, if not miraculous, that a lattice model with 
non-local variables can be described, in the scaling limit, 
by a local field theory. The only trace the lattice non- 
localities leave in the continuum local theory seems to be 
the presence of logarithms in correlation functions. 

It is therefore essential to check that the logarith- 
mic conformal description is indeed appropriate for these 
models, as extensively as it has been done for equilibrium 
critical phenomena (see for instance [7]). 

It is our purpose in this Letter to take further steps 
in this necessary procedure, in the context of the two- 
dimensional Abelian sandpile model. A certain number 
of checks have been carried out for this model (see [4]), 
but the one we propose here is more crucial because it 
deals with microscopic variables which are manisfestly 
non-local, and for which the logarithmic conformal the- 
ory makes a very definite prediction. It therefore exposes 
in the clearest possible way the non-local features of the 
model. 

Namely, we compute, in the infinite discrete plane, the 
2-site correlations Pii{r) — PiPi of two height variables. 



one of which being equal to 1, the other, hi, being equal 
to 2, 3 or 4 (here Pi is the 1-site probability on the in- 
finite plane). Conformal field theory predicts that the 
dominant term of these is given by [5] 

Pu{r) - PiP = c,^^ + . . . , (1) 

with known coefficients Ci . New and explicit lattice calcu- 
lations, to be detailed below, fully confirm these results, 
and exactly reproduces the coefficients Ci. 

Logarithmic conformal theory also predicts that the 2- 
site correlations Pij{r) — PiPj of two heights bigger or 
equal to 2 decay like log^ r/r"*, but the explicit lattice 
calculation of these remains out of range for the moment. 

THE SANDPILE MODEL AND LOGARITHMIC 
CONFORMAL THEORY 

We briefly recall the sandpile model introduced by Bak, 
Tang and Wiesenfeld in [8] (see [9] for further details). 

Every site i of a finite rectangular grid C is assigned 
a height variable hi, taking the four values 1, 2, 3 and 4. 
A configuration C is the set of values {hi} for all sites. 
A discrete stochastic dynamics is defined on the set of 
configurations. If Ct is the configuration at time t, the 
height at a random site i of Ct is incremented by 1, ft.^ — > 
hi + 1, making a new configuration Cj . If the (new) height 
hi in C( is smaller or equal to 4, one simply sets Ct+i = C[. 
If not, all sites j such that their height variables hj exceed 
4 topple, a process by which hj is decreased by 4, and 
the height of all the nearest neighbours of j are increased 
by 1. That is, when the site j topples, the heights are 
updated according to 

hi^hi-Aji, (2) 

with A the discrete Laplacian, An = 4, = —1 for 
nearest neighbour sites, and Ay = otherwise. This 
toppling process stops when all height variables are be- 
tween 1 and 4; the configuration so obtained defines Ct+i- 
The boundary sites are dissipative, because a toppling 
there evacuates one or two grains of sand, which we imag- 
ine are collected in a sink site, connected to all dissipative 
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sites. The presence of dissipativc sites is essential for the 
dynamics to be well-defined, since it makes sure that the 
toppling process stops in a finite time. 

When the dynamics is run over long periods, the sand- 
pile builds up, being subjected to avalanches spanning 
large portions of the system. This correlates the height 
variables over very large distances, and makes the system 
critical in the thermodynamic limit. 

It turns out that, when the dynamics is run for long 
enough, and no matter what the initial configuration is, 
the sandpile enters a stationary regime, in which only 
special configurations occur with equal probability, the 
so-called recurrent configurations [10]. The recurrent set 
TZ forms a small fraction of all configurations, since 

|7^| =detA ~ (3.21)^, (3) 

where N is the number of sites. So the asymptotic state 
of the sandpile is controlled by a unique invariant dis- 
tribution P^, uniform on the set TZ of recurrent configu- 
rations, and zero on the non-recurrent (transient) ones. 
In the infinite volume limit, the invariant measure is 
believed to become a conformal field theoretic measure. 

To be recurrent, the height values of a configuration 
must satisfy certain global conditions [10], leading to non- 
local features. For what follows, it will be enough to 
know that the recurrent configurations are in one-to-one 
correspondence with oriented spanning trees on £*, the 
original lattice C supplemented with the sink site. This 
change of variables, more convenient to perform actual 
calculations, also yields a different lighting on the non- 
localities of the model. 

Spanning trees are acyclic configurations of arrows: at 
each site i of C, there is an outgoing arrow, pointing 
to any one of its A^^ neighbours (if i is dissipativc, the 
arrow can point to the sink site). A configuration of 
arrows defines a spanning tree if it contains no loop. By 
construction, the paths formed by the arrows all lead to 
the sink site which is the root of the tree. 

The mapping between recurrent configurations and 
trees is complicated and non-local; however the spanning 
trees provide an equivalent description. The global con- 
ditions that the heights of recurrent configurations have 
to satisfy are encoded in the property of arrow configura- 
tions of containing no loop, also a global constraint. The 
invariant measure becomes simply a uniform distri- 
bution on the spanning trees. 

Height values at a given site can be related to proper- 
ties of spanning trees. To do so, one defines the notion 
of predecessor: a site j is a predecessor of i if the unique 
path from j to the root passes through i. Then it has 
been shown [11] that the trees in which the site i (not on 
the boundary) has exactly a — 1 predecessors among its 
nearest neighbours correspond to configurations where 
hi > a, for a = 1, 2, 3 or 4. So configurations with hi = 1 
are associated with trees which have a leaf at i: this is 
a local property which may be verified by looking at the 



neighbourhood of i only. In contrast, heights 2, 3 and 4 
correspond to non-local properties in terms of the trees. 

Using this correspondence, joint probabilities for 
heights P[hi = a, hj =&,...] can be related to the frac- 
tions of trees satisfying certain conditions regarding the 
number of predecessors of among their nearest 

neighbours. However, because of the remark we have 
just made, probabilities with heights 1 only are consid- 
erably easier than those involving higher heights. So far, 
the only probabilities involving higher heights in the bulk 
which have been computed are the 1-site probabilities 
P[hi = a] on the upper-half plane [5]. They provided 
enough input to assess the conformal nature of the four 
height variables in the scaling limit. 

The logarithmic conformal theory, relevant to the 
sandpile model, has central charge c = —2. Among the 
distinctive features of a logarithmic theory is the pres- 
ence of reducible yet indecomposable Virasoro represen- 
tations; this property in turn introduces logarithms in 
their correlators [12]. 

The fields describing the scaling limit of the four lattice 
height variables S{hz — i) ~ Pi, which we call hi{z), have 
been determined in [5]. As hinted by the remarks made 
above, the height 1 field is very different from the other 
heights' fields. It turns out that hi is a primary field with 
conformal weights (1,1), while the other three, /i2, ^3 and 
hi, are all related to a single field, identified with the 
logarithmic partner of hi. More precisely, if hi is the 
primary field normalized as the height 1 variable on the 
lattice, then /i2 satisfies the triangular relations, 

Loh2 = h2--hi, Li/?,2 = p, L-ip=--hi, (4) 

where p is a (0,1) field. In fact, hi and /12 are members 
of the non-chiral version of the indecomposable represen- 
tation called 7?,2,i in [13]. The last two fields are linear 
combinations, hs = a3/i2 + Pahi, h^ = 04/12 + (iihi, and 
may also be viewed as logarithmic partners of hi. The 
coefficients ai,Pi are such that /i3 and /i4, like hi and /12, 
have the same normalization as their lattice counterparts; 
their exact values are known [5]. 

The identification of the height fields makes it possible 
to compute correlations. In particular the joint proba- 
bilities for two height variables on the infinite plane 7? 
correspond, in the scaling limit, to 3-point correlators in 
the conformal theory, 

Pij{zi,Z2) - PiPj = {hi{zi)hj{z2)oj{oo)), (5) 

where a; is a weight (0,0) conformal field, logarithmic 
partner of the identity [5]. Indeed the infinite plane 
should be thought of as the limit of a growing finite grid, 
which has dissipation located along the boundary. In the 
infinite volume limit, the boundaries, and with them, the 
dissipation, are sent off to infinity. The field w precisely 
realizes the insertion of dissipation at infinity, required 
for the sandpile model to be well-defined. 
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The 3-point correlators have been computed in [5] , and 
take the general form (212 = 21 — Z2) 



I s, f . . ^, A,j + Bij log Zi2\ + C^j \0g 212 

F12I 

(6) 

where C'ij = if min(z, j) = 1, and moreover Bn = and 
All = —Pil"^ [14], so that, depending on i and one, 
two or three terms in the numerator are present. The 
coefficient of the dominant term, i.e. the largest power 
of log|zi2|, could be determined exactly, and yields the 
dominant contribution of the 2-site probabilities [5] 



Pi,{r)-PiP, 



aiUjPi 



2r4 



log^ r, 



^ > 1, (7) 
«,J>1, (8) 



where Pi = 2(7r — 2)/Tr^, as first computed in [14], and 



^2 = 1, 0-3 = 



8-TT 

2(7r-2 



a4 



2(7r-2)' 



(9) 



CALCULATIONS ON THE LATTICE 



It has been shown in [14] and [11] (see also [5] for 
details) that height probabilities Pi in the ASM can be 
reduced to the computation of determinants of discrete 
Laplacian matrices perturbed by a number of defects. 
The resulting matrices A' = A + B differ from the reg- 
ular Laplacian by a defect matrix B, with B = except 
for a finite number of elements. Given a lattice point to, 
non-zero elements of B related to to can be marked by 
arrows at adjacent bonds (Fig.l). For instance, the non- 
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(c) 




FIG. 1: Non-zero elements of B related to to (arrowed bonds) 
for (a) A' = A + Bi, (b) A' = Aio,,i and (c) A' = Aioop. The 
bond [to,t4] in Aioop is weighted by — s. 

zero part of the matrix B = Bi used for the evaluation 
of Pi (Fig. la) is 



Bi 




(10) 



where rows and columns are labeled by to, t2, ^3, t^. The 
probability to have a height 1 at to is then [14] 

Pi = ^^^^^^^=detil + BiG), (11) 



where G = A ^ . The explicit form of the translation 
invariant Green function on the plane, G(r) = G-^j = 
Go,o+gp,q for f = (p, g). 
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Stt^ JJ_^ 2 — cos a — cos (3 



dad/3, (12) 



implies Pi = 2(7r - 2)/7t^. 

The probability to have a height 2 can be written [11] 

as 



4detAiocai , 4detAioop , 

■^2 ^ Pi-\ \- lim , , . + 

det A £^00 e det A 



+ hm 



2 det A 



e 



[a,6,c] 



£3 det A ' 



(13) 



where the defect matrices related to Aiocai and Aioop are 
shown in Fig.lb-c, and that related to Ae is on the left 
side of Fig. 2. The matrix Ae differs from A by the re- 
moved bond [jo, J3] and three additional matrix elements 
(bonds) weighted by —e between the sites jo, j2, j4 and a 
triplet of neighbouring sites [a, 6, c], whose position and 
orientation (see Fig. 3) are to be summed over the whole 
lattice, with the restriction that the group [a, 6, c] does 
not overlap with jo, j2, is, ji- 

• b 

• c fa 





jo \ 







FIG. 2: Structure of matrix Aie. On the left side, the de- 
fects are the removed bond [ja, jo] and three additional bonds 
[j4,&], [jo,c], [j2,a] with weight —e. The right part shows the 
same defect as in Fig. la. 



c> 



T a. — .c c* — »b '[^ 



FIG. 3: Four possible orientations of the group [a,b,c\. 

The 2-site probability P12 combines the defect of Ai 
with those of Aiocai, Aioop and Ae. Simple calculations 
show that the matrices Aiocai and Aioop contribute a term 
1/r'* to the asymptotics of Pi2{r) for large r, and are 
therefore subdominant. Thus the leading contribution 
comes from the matrix Aie combining the defects of Ai 



4 



and Ae, as shown in Fig. 2. The correlation function 

Pie(r) is 



Pieir) = 

[a,b,c] 



lim 



2 det Aie 



oo £^ det A ' 



(14) 



where the prime means that the sum excludes the terms 
where at least one edge in the group [a, b, c] overlaps a 
deleted edge adjacent to to. The ten forbidden positions 
are shown in Fig. 4. 
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FIG. 4: Forbidden positions of the group [a,b,c]. 

The ratios of determinants in Eq.(14) are computed as 
det(I + BiqG), like Eq.(ll), where the non-zero part of 
BiQ = Aie — A is a 8 X 9 block diagonal matrix. The 
first block is 



-Be = 



1-10 

-e 

-e 

-e 



(15) 



with rows J3, jo, j2, ji and columns jo, js, c, a, 6, while the 
second block is Bi in Eq.(lO). For large r 3> 1, we can 
replace all Green functions containing r by their asymp- 
totic value, 



where we find, after some algebra. 



y9k-i,i-i - 45fc-i,i + gk-i,i+i 

- 9k,i-2 + ^gk,i — gk,i+2 - gk+i,i-i 
+ 'igk+1,1 — gk+i,i+i — ^gk+2,i^- 

The summation over all k, I yields 

+cx> +00 _ 2-j2 



(20) 



(21) 



k= — oo l=—oo 

Finally, we obtain 

P,,(.)_P,P, = _fci^i^ 



+ -r 



(22) 



which coincides with the LCFT prediction (7) for i = 2. 
Similar calculations for Pi^ir) and Pi4{r) fully confirms 
the results (7) with the correct values of the coefficients. 

Despite a very specific form of Ae, the correlation 
functions Pii{r), i = 2,3,4, are the first example where 
the logarithmic corrections to pair correlations can be 
computed explicitly. 
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9.. = ^ ;^U + 4^°s2j' (16) 

where + » 1 and 7 = 0.57721 ... is the Euler con- 
stant. Expanding the determinants in Eq.(14), we ob- 
tain, for the sum over the forbidden positions, 



(17) 



We can now write the desired correlation in the form 



Pie(r) - PiPe = 2 (^Ur{^ - F(r) j , 



(18) 



where the sum is taken over all lattice points s = (fc, I) 
and Pq is the last term in Eq.(13). The function hiris) 



behaves as 



if s > r ^ 1 and as s 



logr, if r » 1, 



s 1, s < r. In the region r 3> 1 and s < r, we have 



(19) 
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